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ABSTRACT 

Theoretical expectations for direct and indirect searches for rehc neutrahnos are 
presented. Complementarity among various investigation means is discussed in 
connection with the values of the neutralino relic abundance. 
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^ ■ 1. Introduction 

^ \ In the present report we discuss the theoretical expectations for detection of rehc 

^ ' neutralinos either by the direct method {e.g. by measuring the nuclear recoil energy 

due to neutralino-nucleus elastic scattering) or by detecting indirect signals due to 
Qh! neutralino-neutralino annihilation occurring in the halo or inside celestial bodies (Sun 

and Earth). 

Evidence for the presence of these relic neutralinos would obviously be of primary 
. importance as a test of the standard Big Bang theory, independently of the fact that 

\ these Susy relic particles would or would not provide a large contribution to the total 

' mass of the Universe. Should the relic neutralinos contribute substantially to the 

total density parameter f2, their detection would also provide an extraordinary hint 
for a clarification of the long-standing dark matter issue. In fact, in this case, relic 
neutralinos would play a significant role as Cold Dark Matter (CDM) constituents. 

In the following we adopt the usual assumptions that: i) R-parity is conserved, ii) 
the neutralino is the lightest supersymmetric particle (LSP) . Very convenient theoret- 
ical frameworks where dark matter neutralino phenomenology may be easily studied 
are provided by the Minimal Supersymmetric Standard Model (MSSM) and by its 
implementation in a Supergravity theory (SUGRA) Here MSSM is meant to 

denote the minimal supersymmetric extension of the standard model where slentons 
and squarks are taken as degenerate, with the exceptions for the stop particles 1^. 

The neutralino (x) is defined as the lowest-mass linear combination of photino, 
zino and higgsinos 
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X = ai^ + a2Z + a^iHl + 04^^° (1) 

Here 7 and Z are the fields obtained from the original U(l) and SU(2) neutral gaug- 
inos, B and by a rotation in terms of the Weinberg angle. 

The neutralino mass and the coefficients Oj depend on the parameters: /x (Higgs 
mixing parameter), Mi, M2 (masses of B and of W3, respectively) and tan/3 = Vu/vd 
{vu and Vd are the v.e.v.'s which give masses to up-type and down-type quarks). It 
is customary to employ the standard GUT relationship between Mi and M2: Mi = 
(5/3) tan^ 6wM2 — 0.5 M2. We use this assumption here. 

In the following for the parameters M2 and /i we will consider the ranges: 
20 GeV< M2 < 6 TeV, 20 GeV< \n\ < 3 TeV. tan/? will be taken at the representa- 
tive value tan (3 = 8. 

For the evaluation of the neutralino relic abundance and of the event rates for 
direct and indirect neutralino detections one has also to assign values to the masses 
of a large number of particles, namely to the Higgs bosons and to the Susy scalar 
partners of leptons and quarks: sleptons (/) and squarks (g). In the MSSM scheme 
we consider here these values are assigned arbitrarily: a standard procedure consists 
in assuming mass degeneracy both for sleptons and for squarks except for the stop 
particles As for the neutral Higgs bosons we recall that in the MSSM there are 
three neutral Higgs particles: two CP-even bosons h and H (of masses rrih, itlh 
with rriH > mh) and a CP-odd one A (of mass m^). Once a value for one of these 
masses (say, rrih) is assigned, the other two masses (m^, uih) are derived through 
mass relationships depending on radiative effects. 

Implementation of MSSM with supergravity sets a much more constrained phe- 
nomenological framework, since SUGRA establishes strict relations between all the 
masses in play and the few fundamental theoretical parameters: A and B (appearing 
in the soft symmetry-breaking interaction terms), ttlq (common scalar mass at the 
GUT scale), mi/2 (common gaugino mass at the GUT scale) and /i. Furthermore, 
other specific theoretical requirements (features of the symmetry breaking, condition 
that the neutralino be the LSP, ...) strongly restrict the whole parameter space. This 
has important consequences; for instance, it constrains the neutralino to composi- 
tions with dominance of the gaugino components. It has to be noted that important 
constraints on the neutralino parameter space may be inferred from the recent data 
on 6 — > s + 7 process §. The precise nature of these constraints still requires further 
investigation, especially in view of the large uncertainties due to QCD-effects I. 

In the present note the emphasis is on the investigation power of various detection 
methods and on their complementarity. On purpose, the model adopted here for 
the calculations is as simple as possible (not too much constrained by theoretical 
requirements) with choices for the free masses that are only restricted by experimental 
bounds. 

Our main concern is to discuss the minimal sensitivity required in experimental 
devices in order to undertake a significant investigation of neutralino dark matter. 
For this reason we present here evaluations where the smallest values compatible with 



experimental lower bounds are assigned to the unknown masses; this usually provides 
maximal values for the signals. To be definite, in the following we will set the sfermion 
masses at the value = 1.2 m^, when > 45 GeV, = 45 GeV otherwise. Only 
the mass of the top squarks are assigned a larger value of 1 TeV. The Higgs mass 
is set at the value of 50 GeV. The top mass is fixed at nit = 170 GeV. 



2. Neutralino Relic Abundance 



For the computation of the direct and indirect event rates for neutralino one has to 
use a specific value for the neutralino density p^. Obviously, it would be inappropriate 
to assign to the neutralino local density the standard value for the total dark matter 
density pi = 0.3 GeV cm~^, unless one specifically verifies that the neutralino relic 
abundance turns out to be at the level of an (f2/i^)min consistent with pi. This is 
why a correct evaluation of the event rates for x detection also requires a calculation 
of its relic abundance. 

Thus we evaluate f^^^^ determine p^ by adopting a standard procedure o-. 

when Q^h"^ > {Qh'^)^i^, we put p^ = pi; when Q^h"^ turns out to be less than {Qh'^)^i^, 
we take 

^""^^ (2) 

(Vlh )min 

Here (f2/i^)mm is set equal to 0.03. 

For the neutralino relic abundance fi^^^ employ the results of our previous 
work i. In Fig.l we display regions of the M2,p plane which are characterized by 
different values of Q^h"^. Also shown in this figure are the iso-mass curves (dashed 
lines) and the iso-compositions curves (solid lines). Along an iso-composition line the 
composition parameter P defined as the gaugino fractional weight, i.e. P = + a^, 
is kept fixed. In this figure, as well as in the following ones, only results for positive 
p are displayed. 

The results of Fig.l can conveniently be reported in a ^^^"^ "^x P^ot, for 
fixed values of P. This is done in Fig. 2 where fi^^/i^ is plotted as a function of 
for three representative neutralino compositions: i) a gaugino-dominated composi- 
tion (P = 0.9), ii) a composition of maximal gaugino-higgsino mixing (P = 0.5), 
iii) a higgsino-dominated composition (P = 0.1). As expected, out of the three 
compositions displayed in Fig.2, the gaugino-dominated state provides the largest 
values of Q^h"^. In order to have more substantial values of Q^h"^, one has to con- 
sider purer gaugino compositions (P > 0.99). The very pronounced dips in the plots 
of Fig.2 are due to the s-poles in the x^X annihilation cross section (exchange of 
the Z and of the neutral Higgses). The sharp decrease at 80-90 GeV is due to the 
opening of the W~^W~ and ZZ final states in the x^X annihilation. We remind that 
^x^"^ ^ '^a^nnV >int)~\ where < cTann'f >int IS the annihilation cross section times 
relative velocity, averaged over the neutralino thermal distribution, integrated from 
the freezing temperature down to the present temperature. 



3. Indirect Detection at Neutrino Telescopes 



Let us turn now to the indirect search for neutralino dark matter which can be 
performed by means of neutrino telescopes Neutralinos, if present in our 

galactic halo as dark matter components, would be slowed down by elastic scattering 
off the nuclei of the celestial bodies (Sun and Earth) and then gravitationally trapped 
inside them. Due to the process of neutralino capture these macroscopic bodies could 
accumulate neutralinos which would subsequently annihilate in pairs. An important 
outcome of this x^X annihilation would be a steady flux of neutrinos from these 
celestial bodies. 

The differential neutrino flux at a distance d from the annihilation region is given 

by 



dE„ ~ And^ h dE, 



E^i/'^ (3) 



where T a is the annihilation rate and F denotes the x^X annihilation final states 
which are: 1) fermion-antifermion pairs , 2) pairs of neutral and charged Higgs bosons, 
3) one gauge boson-one Higgs boson pairs, 4) pairs of gauge bosons; -B^j"* denotes 
the branching ratio into the fermion / (heavy quark or r lepton), in the channel 
F] dNfy/dEy denotes the differential distribution of the neutrinos generated by the 
semileptonic decays of the fermion /. The i/^'s, crossing the Earth, would convert 
into muons and generate a signal of up-going muons inside a neutrino telescope. 
Calculations of this muon flux from the original neutrino flux may be performed 
using standard procedures 

Particular care has to be taken in the evaluation of the annihilation rate Ta- This 
quantity is given byi 

TA = ^t^n\i^{—] (4) 



2 \ta 

where t is the age of the macroscopic body (t = 4.5 Gyr for Sun and Earth), ta = 
(CC^)~^/^, C is the capture rate of neutralinos in the macroscopic body and Ca is 
the annihilationrate per effective volume of the body. The capture rate C is provided 
by the formula EJ 

C=P^Y.^(Msh){vl^hX,^ (5) 

where is the neutralino mean velocity, a^^i is the cross section of the neutralino 
elastic scattering off the nucleus i of mass rrii (for some properties of the elastic 
nucleus cross section see next Sect. 4 and Ref. 15), Msfi is the total mass of the 
element i in the body of mass M^, {vl^^i is the square escape velocity averaged over 
the distribution of the element Xj is a factor which takes account of kinematical 
properties occurring in the neutralino-nucleus interactions. Ca is given by 



where < aannV >o is the annihilation cross section times relative velocity, averaged 
over the neutralino thermal distribution, at present temperature. Vq is defined as 
Vo = {3mliT/{2p X 10 GeV))^/^ where T and p are the central temperature and the 
central density of the celestial body. For the Earth (T = 6000 K, p = 13 g ■ cm~^) 
Vo = 2.3xl025cm^ for the Sun (T = 1.4xlO^K, p = ISOg-cm"^) ^ = 6.6x10^8 cm^. 

For the computation of the capture rate (and then also of r^) one has to use a 
specific value for the neutralino density p^ according to the procedure explained in 
Sect. 2. 

The results of our evaluations for the quantities C, Ca and are reported in 
Ref. 11,12. From these results it turns out that: i) for the Sun, the equilibrium 
between capture and annihilation is reached over the whole range; ii) in the case 
of the Earth, equilibrium is not reached for > mw, because of the substantial 
suppression introduced in F^ by the factor tanh^(t/r^). 

Now we report some of our results about the fiux of the up-going muons in the 
case of x^X annihilation in the Earth. In Fig. 3 we show the fiuxes of the up-going 
muons as functions of for a number of values of the neutralino composition P, for 
X^X annihilation in the Earth. The threshold for the muon energy is E^^ = 2 GeV. 
We recall here that the present experimental upper bound for signals coming from 
the Earth is 4.0 • 10~^^cm~^s~^ (90 % CL.) 0. By comparing this upper limit with 
our fiuxes we see that the regions explored by Kamiokande (at our representative 
point: tan/3 = 8,mh = 50 GeV) concern the mass range 50 GeV < < 65 GeV. 
These regions are illustrated in Fig. 4 in a M2-/i plot. In this figure we also display 
the regions which could be explored by a neutrino telescope with an improvement 
factor of 10 (and of 100) in sensitivity. 

The location and the shape of the most easily explorable regions in the M2-/i plane 
depend on the Earth chemical composition and on the neutralino composition in terms 
of the gaugino, higgsino components. In fact the capture rate of neutralinos is more 
effective when neutralino mass matches the mass of some of the main components 
of the Earth and when the neutralino is a large gaugino-higgsino mixture. Because 
of these two properties the signal is maximal along iso-mass lines in the range 50- 
65 GeV, with elongations along iso-composition lines of sizeable mixing. 

In Fig. 5 we report the fluxes for up-going muons as functions of due to x^ 
X annihilation in the Sun. As before the threshold for the muon energy is E^^ = 
2 GeV. The evaluated fluxes are below the present experimental upper limit of 
Kamiokande: 6.6 ■ 10~^^cm~^s~^ (90 % C.L.) 11^. The regions explorable by a neutrino 
telescope with an improvement factor of 10 (and of 100) in sensitivity respect to 
the present Kamiokande sensitivity are displayed in Fig.6. Here the regions that 
are more easely explorable, expand toward gaugino-dominated compositions since 
spin-dependent cross sections (with exchange of light squarks) are important in the 
capture of neutralinos by the Sun. 



From these results it can be concluded that neutrino telescopes with an area 
above 10^ are very powerful tools for investigating neutralino dark matter in large 
regions of the parameter space. It also emerges from the previous results that the 
signals from the Earth and from the Sun somewhat complement each other to allow an 
exploration about DM neutralino over a wide range of m^. In Ref. 12 has been derived 
the relation between the exposure At of a neutrino telescope {A being the telescope 
area, t the live time) and the explorable range in the neutralino parameter space, 
when the signal-to-background ratio is optimized by appropriate angular selections. 



4. Direct Detection 



Another way to search for dark matter neutralinos is the direct detection which 
relies on the measurement of the recoil energy of nuclei of a detector, due to elastic 
scattering of x's. The relevant quantities to calculate are the differential rate (in the 
nuclear recoil energy Er)- 

-Nt^ dvf{v)v—^{v,Er) (7) 



dEj. Jvmin{Er) dEj. 

and the integrated rate -Rint, which is the integral of Eq. (7) from the threshold energy 
E^^, which is a characteristic feature of the detector, up to a maximal energy E^^^. 
In Eq.(7) denotes the number of target nuclei, daei/dEr is the differential elastic 
cross section and f{v) is the distribution of x velocities in the Galaxy. It is important 
to note again that the local density is evaluated here according to the procedure 
discussed in Sect. 2. In general, the x^nucleus cross section has two contributions: a 
coherent contribution, depending on A"^ {A is the mass number of the nucleus) which 
is due to Higgs and q exchange diagrams; a spin-dependent contribution, arising from 
Z and q exchange, proportional to }?J{J + 1). 

By way of example, let us remind the expression of the coherent cross section due 
to the Higgs-exchange El: 

CelH — «H-A 



Hi A f I \o \ ' 

TT [rrii + m^y 

where an is a quantity depending on the neutrahno-Higgs and the Higgs-quarks 
couplings. It is worth mentioning that an depends rather sensitively on the 
composition and on a number of parameters, such as tan (3 and the Higgs masses. 

Except for very special points in the parameter space, the coherent contribution to 
elastic cross section strongly dominates over the spin-dependent one. For a detailed 
analysis on the calculation of the direct event rates see Ref. 19 and references quoted 
therein. For an experimental overview about dark matter detectors see Ref. 20. 

Here, as an example, we simply report in Fig. 7 the event rates R\-at for a Germa- 
nium detector as a function of for neutralino compositions P = 0.1, 0.5, 0.9. Rates 
are calculated by integrating the differential rate of Eq.(7) over the electron-equivalent 
energy range (2-4) KeV. In Fig. 8 we show the regions of the M2-/i parameter space 



which can be explored with an improvement of one and two orders of magnitude in 
the sensitivity of the detector. 

As for the shape of these regions we refer to the comments presented above, in 
Sect. 3, in connection with Fig. 4. Again, the signals are higher along the iso-mass 
line with an close to the mass of the nuclei composing the detector. Thus, using 
detectors of different compositions allows explorations of the M2-/i plane over a wide 
range in m~^. For instance, investigation of regions with small rriy. values (of order of 
10 GeV) with very low threshold detectors c3 would be very interesting. In fact this 
rriy. range (which is excluded by accelerator data only under a number of assumptions) 
is out of reach for the indirect detection discussed in the previous Section. 

5. Comparison between Direct and Indirect Signals 

Fig.s 3,5,7 provide a comparison between the capabilities of the direct detection 
versus the detection by neutrino telescopes. It has to be noted that these figures refer 
to an arbitrary representative point of the parameter space. However, since both 
(direct and indirect) signals depend mainly on a common quantity, i.e. the x^nucleus 
elastic scattering, the relative size is approximately independent of the variation of 
the model parameters. The main departures from this feature occur in the case of 
the signal from the Earth in regions where the factor tanh."^ (t / ta) is small and for 
the signal from the Sun in regions where spin-dependent effects in the capture rate 
are dominant. Fig.s 3,5,7 also display how the the theoretical predictions compare 
with the present sensitivity for each individual method. A comparison of the relative 
power of the two detection methods for exploring the M2-/X plane is provided by Fig.s 
4,6,8. 

6. Conclusions 

In the previous sections we have examined two types of detection methods for relic 
neutralinos, one based on direct measurement and one on the detection of indirect 
signals due to x^X annihilation in celestial bodies. As we have seen by the explicit 
expressions given in Sect.s 3-4, in both cases the detection rates are proportional 
to the neutralino local density and to the elastic neutralino-nucleus cross section: 
R oc Px^ei (for simplicity, here and in the following, we consider the neutralino mass 
fixed at some arbitrary value m^; the other free parameters may vary in the parameter 
space along iso-mass curves). 

Because of the properties discussed in Sect. 2, once the neutralino relic abundance 
is evaluated, one has to distinguish between two cases: 

a) — ^min which entails R oc (Xei 

b) fly. < flmin which entails R oc ael/ < CTannV >int 

Now, if the strength of the couplings increases, and hence aei and (Tann increase, 
it turns out that also in case b) the rates R usually increase even if at a much less 
extent than in case a). It then happens that for the two detection methods discussed 



so far, the largest values of the detection rates for relic neutralinos occur for models 
where tan/3 is large and neutralino is largely mixed (this last condition has to be 
somewhat relaxed for the signal from the Sun, as discussed above). For these models, 
because of the large value of the annihilation cross section, the neutralino relic density 
is small. This is a scenario where the detection methods under discussion have good 
chances to detect relic neutralinos, even if these particles cannot provide a substantial 
contribution to the total Q of the Universe. Regions of the parameter space where 
this situation occurs are defined by some authors as cosmologically uninteresting and 
consequently disregarded. On the contrary, we believe that they deserve a careful 
investigation, since much interesting physics could occur there. The relevance of this 
point is further confirmed by some recent results derived from SUGRA theories where 
some of the universality properties, usually required at the GUT scale, are relaxed B3. 
In fact it turns out that, in this case, very interesting schemes with large neutralino 
mixing, large tan (3 and small relic abundance naturally emerge. 

As we have seen, when is large the two detection techniques previously dis- 
cussed have less chances to succeed (except for the signal from the Sun, which however 
requires neutrino telescopes of very large area). In order to cover this case adequately, 
one may resort to different detection methods. These are based on the measurements 
of indirect signals due to neutralino-neutralino annihilation in the Galactic halo. 
These have been widely discussed in the literature; we refer to Ref. 23 for previous 
references and for a new analysis aimed at the evaluation of the antiproton to proton 
ratio (p/p) in cosmic rays. A common feature of this class of measurements is that 
the event rates are typically proportional to the square of the local neutralino density 
times the annihilation cross section: R oc < CTann^ >o- This implies that, in this 
case, detection rates are high in regions of the parameter space where the neutralino 
relic abundance is large. By way of example, we give in_Fig.9 the regions of M2-/i 
plane that can be explored by measuring the (p/p) ratio c3. 

In conclusion we may state that complementarity among various detection meth- 
ods may potentially offer a good coverage of the neutralino parameter space. 
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